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in our study. The first one concerns the transient process of self-induced transparency 
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collisions of counter-propagating solitons. It is shown that the level of inhomogeneous 
broadening has substantial effect on elasticity of such collisions. 
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Self-induced transparency (SIT) is one of the basic phenomena of nonlinear optics. This 


effect resulting in formation of temporal solitons (so-called 27r pulses) in a resonantly 


absorbing medium was described in the pioneering papers by McCall and Hahn m- 


There is an extensive literature on this topic including the classic monograph by Allen 


and Eberly [3] and a number of reviews (see @1 El El El, to name a few). The continuing 
study of SIT is motivated by the fundamental and general character of the two-level 
model which is the basic tool for explanation of this phenomenon. The investigation 
of the properties of this model based on the semiclassical Maxwell-Bloch equations 
is necessary for our understanding of the nature of light-matter interaction. There 
is a number of generalizations, as well, aiming to take into account the near-dipole- 
dipole interactions (local-held correction) in the dense resonant media |8], Stark shift of 
the absorption line [9], and few-cycle pulse dynamics in the regime of invalidity of the 
rotating-wave approximation mm • Using the two-level model and its generalizations, 
the broad spectrum of SIT studies was performed including invariant pulse propagation 
and optical switching in dense media [121 [13], quasiadiabatic following analysis [14] , 
incoherent solitons [15], SIT soliton lasers [16], SIT soliton collisions [T7, HE], SIT in 
Bragg rehectors and photonic crystals [19) [20, [21], coherent pulse propagation and 
SIT effects in doped waveguides and amplifiers [221, 23], SIT in the presence of Kerr 
nonlinearity [2D [25], etc. 

The aim of this paper is to consider some details of SIT soliton formation and 
interaction between solitons in the two-level medium with inhomogeneous broadening 
of resonant line. This broadening is due to the fact that, generally, the frequency 
of resonant transition is not the same for all atoms of the medium. The importance 
of inhomogeneous broadening can be illustrated by the example of the so-called area 
theorem. This theorem governs the change of the “area” S(z ) = (2 fi/h) f A(z,t)dt of 
the pulse propagating in the medium and can be written in the form 



( 1 ) 


where /i is the component of transition dipole moment parallel to the polarization vector 
of the electric held, A the electric held amplitude, So the initial area, a the absorption 
coefficient of the medium, h the Planck constant. It is known that the area theorem 
(Jl]) is strictly valid only for the inhomogeneously broadened media, though the main 
features of SIT can be observed in the case of homogeneous broadening as well [6]. This 
was confirmed recently by Yu et al. [26] who performed the direct simulations of area 
change. They demonstrated that in the case of homogeneous broadening, approach of 
the area to the stationary value corresponding to the SIT soliton is accompanied by 
slowly damping oscillations. These oscillations cannot be described by the standard 
formulae such as 

The present paper continues previous studies of inhomogeneously broadened media 
in comparison to their homogeneously broadened counterparts. Our attention is directed 
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to the two questions studied in our recent works [TH], [27J in the approximation of 
homogeneous broadening, viz. the transient process of SIT soliton formation and the 
collisions of such solitons. After Sec. [2] where the model is described, we discuss these 
two questions in Sec. [3] and 0 respectively. The paper closes with the brief conclusions 
of the results obtained. 


2. Problem statement 


Light interaction with the two-level medium at every spatial point is governed by the 
system of semiclassical Bloch equations for population difference W and microscopic 
polarization R : 

HR 

— = iSlW + iR6 - 72 i?, (2) 

dr 

dW 

— = 2i(n*R- R*n) » 7 i(IP-l), (3) 

dr 

where hi = (/ i/hu)A is the normalized electric held amplitude (or dimensionless Rabi 
frequency), and 71 = (u/Ti ) -1 and 72 = (uT 2 )^ 1 are the rates of longitudinal and 
transverse relaxation, respectively. Since we are interested in consideration of the 
inhomogeneously broadened medium, the variables W and R are the functions of the 
normalized detuning S = 5 a + Au/u, which is the sum of S a = (u—u a )/u, the normalized 
detuning of the held carrier (central) frequency c a from the average atomic resonance, 
and the term Auju which describes the deviations of the resonance frequency from the 
average value. The distribution of two-level atoms over the detunings is governed by 
the weight function g(5), so that the amplitude of the macroscopic polarization of the 
medium is given by 


p = pC 


R(S)g(5)dS, 


(4) 


where C is the concentration of the two-level atoms. 

The one-dimensional Maxwell wave equation for light pulse propagation in 
resonantly absorbing medium is 

8 2 E 1 d 2 E _ 4tt d 2 P 
dz 2 c 2 dt 2 c 2 dt 2 ’ 

Assuming E = A exp [i(ut — kz)\ and P = pexp[i(ut — kz )], this equation can be 
represented as the expression for the dimensionless held amplitude, as follows [18] 


d 2 n 

W 


d 2 n on on 

— + 21 — + 2 '— 


= 3e 


d 2 p dp 

8t 2 2l dr P 


( 6 ) 


where e = cul/cu = Anp 2 C is the normalized Lorentz frequency which describes the 
strength of light-matter coupling, and p is the integral in the right-hand side of ( 0 ]). 
In all the equations above, t — ut and £ = kz are dimensionless time and distance, 
respectively; k = u/c the wavenumber; and c the light speed in vacuum. Here we 
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assume, without loss of generality, that the background dielectric permittivity of the 
medium is unity, i.e. we consider the two-level atoms in vacuum. It is also important 
to note that in (jBJ) the slowly varying envelope approximation (SVEA) is not used. 

To solve ©-© self-consistently, we apply the numerical approach which is 
essentially the same as in our previous studies EH- At the edges of the calculation region, 
we apply the so-called absorbing boundary conditions using the total held / scattered 
held (TF/SF) and the perfectly matched layer (PML) methods [28j[29]]. One should then 
supplement the method, at every time step, with the solution of the Bloch equations for 
different detunings with subsequent calculation of the integral p according to (141). Note 
that, according to the numerical method, we deal with the evolution of the total electric 
held (or complex amplitude which contains all the changes of phase due to propagation), 
the rehected and transmitted helds being calculated at the corresponding boundaries 
due to the TF/SF procedure. The incident held and the direction of propagation are 
initialized at these boundaries as well. In other words, we do not divide the held into 
two counterpropagating waves inside the medium, in contrast to the frequently used 
approach (see, for example, [30i 3T1132]). 

It is known that in the SVEA regime, there are analytical solutions of coherent 
pulse propagation obtained with the inverse scattering theory (33, ,3U [35]. In particular, 
this method allows to derive the stationary pulse prohle, the temporal dynamics of 
pulse area and effect of relaxation. However, the numerical simulations cover broader 
spectrum of situations, e.g. they allow to take into account the local-held correction 
or deviations from slowly varying envelope and rotating-wave approximations. The 
same can be said about collisions of SIT solitons: as far as we know, there are no 
analytical solutions of this problem. Therefore, in this paper, we use only the numerical 
approach and leave the analytical attempts for future considerations. Our choice of 
one-dimensional approximation (rather than full 3D simulations) is justihed by the 
possibility of controlling instabilities and other effects of transversal beam structure 
by the manipulations with the aperture of the optical system [36]. 

Let us discuss the main parameters of calculations. To preserve the generality, 
all the values are represented in dimensionless form, the central wavelength of pulse 
A = 2hc/uj being the main parameter of normalization. We deal with the ultrashort 
pulses of Gaussian shape D = Q p exp(— t 2 /2tp), where t p is the pulse duration; such 
expressions are used as boundary conditions at the left (z = 0) and right (z = L) 
boundaries of the medium according to the TF/SF method. Throughout the paper, the 
amplitude of the pulses is measured in the units of the characteristic Rabi frequency 
Do = A/2\/27r ct p , which corresponds to the area equal to 2n. The pulse duration is 
measured as a number N of periods of electric held oscillations T = X/c giving its full 
width at half maximum (FWHM), namely t p = NT /2\J\n 2. Since we are interested 
in the study of coherent interaction of light with the resonant medium, we assume 

7 1 = 72 = 0 in our calculations, i.e. the homogeneous broadening governed by qx and 

72 is absent. 

For the distribution function describing the inhomogeneous broadening, we use the 
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Figure 1. Intensity profiles of the 37 t pulse transmitted through the two-level medium 
with different levels of inhomogeneous broadening. The thickness of the medium 


L = 100A. 


Gaussian envelope 

9 ( 8 ) = -i-exp(-<5 2 /<5jj), ( 7 ) 

y/ndo 

where 5 0 = Auio/o; is the normalized width of the distribution, the value Aa; 0 is inversely 
proportional to the characteristic relaxation time T 2 . Such distributions as that in (J7|l 
appear, for example, as a result of Doppler broadening in atomic gases. It is convenient 
to express the frequency parameters through the pulse duration as it was done above for 
the amplitude O 0 - Therefore we adopt 5 a = 0 (exact resonance in homogeneous case) 
and uiLtp = 0.1. The latter condition also allows us not to take into account the so- 
called local-field correction [8] which can be neglected, when Do A (u L t p ) 1 > 1 [37]. 
Similar normalization can be done for the width of inhomogeneous broadening which is 
defined through the value A 0 = Aco 0 t p , so that the condition A 0 <1 (or, equivalently, 
t p -C T 2 *) corresponds to the case of the homogeneously broadened two-level medium. 
The concrete values of parameters can be obtained if one assumes specific values of the 
wavelength A and the number of cycles N. 

3. Soliton formation 

First of all, let us consider the influence of inhomogeneous broadening on the dynamics of 
self-induced transparency (SIT) soliton formation. We take very short pulses containing 
only N = 25 cycles. Such small duration is convenient from the standpoint of 
calculation speed, while, at the same time, it is not too short to break the rotating-wave 
approximation used in the Bloch equations. Anyway, the results can be rescaled for 
other values of t p . As an example, let us consider the propagation characteristics of the 
37 t pulse (with amplitude Q p = 1.5D 0 )- Figured] shows the results of transmission of such 
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Figure 2. (Colour online) Peak intensity of the 37r pulse as a function of distance 
traversed through the two-level medium with different levels of inhomogeneous 
broadening. 


a pulse through the two-level medium of thickness L = 100A. This figure demonstrates 
some important aspects to be mentioned here: (i) the influence of inhomogeneous 
broadening appears when A 0 ~ 1, (ii) increasing A 0 leads to the rise in pulse speed 
in conformity with the previous studies |26j . These facts imply that our method of 
calculation works well, so we can directly proceed to the topic of this paper. 

The dynamics of soliton formation are traced in Fig. [2] where the change in peak 
intensity of the pulse as it propagates in the medium is depicted. As previously, the 
incident pulse has the area 37T. The final, quasistationary state (the intensity of the 
formed soliton) is almost identical in homogeneous and inhomogeneous broadening 
cases, while the initial behaviour differs strongly. The pulse in homogeneously broadened 
medium (A 0 = 0) experiences strong compression soon after incidence and then, through 
sharp oscillations of peak intensity, reaches quasistationary, solitonic form. These 
oscillations are analogous to those in the area dynamics studied by Yu et al. [26] 
and are characteristic for this case. Introduction of inhomogeneous broadening (see the 
curve for A 0 = 1) makes these oscillations less pronounced, while, for A 0 1, they are 
entirely smoothed out. One can conclude that for smaller broadenings, the pulse can be 
stronger compressed during the transient process of soliton formation. 

As another measure of inhomogeneous broadening, we will use the mean value of 
population difference calculated the same way as p, namely 

W = j W(S)g(S)d5. (8) 

In the limit of homogeneous broadening, when all the atoms have the same resonant 
frequency, this integrated population difference gives simply population difference IF. 
Under the influence of, say, 2tt pulse, IF demonstrates typical cycle of inversion depicted 
in Fig. |3](a) : beginning at the ground state (IF = 1), the medium switches to the fully 
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p 

Figure 3. Dynamics of the integrated population difference W at the entrance of 
the two-level medium with different levels of inhomogeneous broadening. The incident 
pulse has the area 2n. 

inverted state [W = —1) and subsequently returns to the ground one (W = 1). The 
same behaviour is seen in Fig. [3(b) for A 0 = 0.1 which is effectively still homogeneous 
case. As the inhomogeneous width of the spectral line grows, the energy of the pulse 
distributes over atoms with different resonant frequencies, so that the integral value 
W cannot reach the full inversion as shown in Fig. [3(c) for A 0 = 1. For strong 
inhomogeneous broadening [A 0 1, see Fig. [3(d)], W stays near the ground state 

at every time instant. In the next section, we will discuss the implications of these 
population difference dynamics for collisions of solitons. 

4. Collisions 

In this section, we study influence of inhomogeneous broadening on the inelastic 
collisions of counter-propagating pulses in two-level medium. We are especially 
interested in the so-called asymmetric collisions when two colliding solitons are not 
identical. It was shown previously [18] that in the case A 0 = 0, such a collision can 
result in total destruction of one of the solitons if the initial parameters of the pulses 
are chosen properly. Since this effect is accompanied by strong absorption of light in 
the point of collision, it was called the controlled absorption of the soliton and used as 
a source of diode action m But what if Ao is not zero? 

Let us consider the collision of two pulses propagating in the two-level medium 
of thickness L = 100A: the first, conventionally called forward propagating (FP, from 
left to right), is the 2tt pulse = fl 0 ), while the second, backward propagating (BP, 
from right to left), is the 37T pulse = 1.5fi 0 )- The results of collisions (profiles of 
FP and BP transmitted radiation) for different levels of inhomogeneous broadening are 
shown in Fig. [4] One can see that at low broadening's (A 0 = 0.1 and 1), the inelastic 
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Figure 4. Intensity profiles of (a) the 2n forward propagating (FP) and (b) the 37 t 
backward propagating (BP) pulses after the collision inside the two-level medium with 
different levels of inhomogeneous broadening. The thickness of the medium L = 100A. 
The profiles are recorded at right and left boundaries of the medium, respectively. 


collision results in total breakdown of the FP 2ir pulse: there is no soliton at the exit 
of the medium. The radiation contains low-intensity oscillations, while most part of 
energy is trapped by the medium around the point of collision as seen in Fig. [5] where 
the distribution of the integrated population difference W after the collision is plotted. 
For larger values of Ao, the FP soliton appears at the exit and its intensity grows with 
increasing A 0 as is clearly seen in Fig. |4])a). Simultaneously, the excitation of medium 
diminishes: obviously, the strong inhomogeneous broadening cannot provide interpulse 
interaction large enough to effectively trap radiation. It is also should be noted that the 
BP 37 t pulse also loses part of its energy due to collision, since the peak intensity drops 
from about 4.5f2o ( see Fig. |T]or m to approximately 30§ (at A 0 = 0.1) and 2.50§ (at 
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Figure 5. Distribution of the integrated population difference W along the two-level 
medium after the collisions depicted in Fig. [I] The distributions are plotted at the 
instant 100i p . 

A 0 = 3) and grows again for larger broadenings (see Fig. 01(b)). Nevertheless, the BP 
pulse is always present at the exit of the medium and can be considered as means for 
controlling less intensive FP soliton. 

These relationships can be traced in more detail in Fig. El where one can see the 
dependencies of the part of the FP pulse energy transmitted through the medium and 
of the peak intensity of the transmitted soliton on the parameter of inhomogeneous 
broadening A 0 . It should be said that the transmitted energy is present not only in the 
form of soliton but also as unstructured radiation (dispersive waves). The difference 
can be clearly seen since the SIT solitons always have characteristic shape described by 
hyperbolic secant. The curves for the two variants of light-matter coupling are shown 
in Fig. El First, let us discuss the case of strong coupling ( u>Lt p = 0.1) corresponding 
to the results represented in Figs. [4] and [5] The soliton is absent at the exit at low 
values of broadening (A 0 < 2.5), though the substantial part of its energy (up to 50%) 
is transmitted in the form of low-intensity oscillations like those discussed above. The 
minimal value of transmitted energy (only about 6%) is observed at Ao — 2.2. Further 
increase in broadening leads to the rapid jump of transmitted energy corresponding to 
appearance of the single solitary pulse at the exit. Subsequently monotonous growth 
of transmitted energy and peak intensity of soliton occurs as is clearly demonstrated in 
Fig. El 

Weakening the light-matter interaction (ajptp = 0.01) results in the pronounced 
shrinkage of the range of broadenings where the transmitted soliton is absent. At the 
same time, the minimum of transmitted energy shifts to the lower value of Ao and 
has significantly larger magnitude. At large broadenings, the pulses experienced the 
collision are more intensive than in the previous case and have the parameters close 
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Figure 6. Dependence of (a) part of energy of 2n pulse transmitted in forward 
direction (by the instant 100i p ) and (b) the peak intensity of transmitted soliton on 
the level of inhomogeneous broadening Ao. The curves are shown for two strengths of 
light-matter coupling. 



Figure 7. The same as in Fig. [61(a) but for 37r pulse transmitted in backward direction. 


to the nonperturbed solitons. Finally, Fig. [7] demonstrates similar dependencies of the 
transmitted energy for the BP 37T pulse. It is seen that increase in the light-matter 
coupling results in reinforcement of inelasticity, since the pulse loses considerably more 
energy due to collision at uiLtp = 0.1 than in the case of UL,t p = 0.01. It should be also 
noted that the position of the minimum at c 0 ]J, p = 0.1 approximately coincides with 
that of Fig. Eta) ; it is not so evident for = 0.01. Nevertheless, one can say that 
there exists the optimal inhomogeneous broadening for observation of inelastic soliton 
collisions. 
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This study of coherent pulse propagation in the two-level medium is focused on the 
influence of inhomogeneous broadening on the SIT soliton formation and on the collisions 
of solitons. Though the main features of these processes remain unchanged, there 
are some interesting details which distinguish the case of strongly inhomogeneous 
broadening from its homogeneous counterpart. Formation of the SIT soliton in the 
homogeneously broadened medium is accompanied by the characteristic oscillations of 
pulse intensity (and its area), so that there is a distance of optimal compression of the 
pulse. At this distance, the pulse is more intensive than after ending of the formation 
process. In strongly inhomogeneously broadened case, these oscillations are entirely 
washed out and, hence, the pulse cannot be compressed more than it occurs in the final 
state of SIT soliton. The inelastic collisions of counter-propagating SIT solitons has the 
optimal conditions as well: there is a nonzero level of inhomogeneous broadening which 
provides the optimal trapping of radiation by the medium due to the collision. Further 
increasing the width of broadened line or decreasing the light-matter coupling strength 
(e.g. due to diluting the medium) results in sharp rising of elasticity of collisions. We 
believe that the results of this study will be helpful for understanding of some subtle 
details of light-matter interaction. 
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